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ABSTRACT 

We present comprehensive cluster membership and g'r' photometry of the 
prototypical old, metal-rich Galactic star cluster NGC 6791. The proper-motion 
catalog contains 58,901 objects down to g' ~24, limited to a circular area of radius 
30'. The highest precision of the proper motions is 0.08 mas yr^^. Our proper 
motions confirm cluster membership of all main and also some rare constituents 
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of NGC 6791. The total number of probable cluster members down to g' = 22 
(My ~ +8) is ~4800, corresponding to Mtot ~ 5000 Mq. New findings include 
an extended horizontal branch in this cluster. The angular radius of NGC 6791 
is at least 15' (the effective radius is Rh — 4.4' while the tidal radius is ~ 23'). 
The luminosity function of the cluster peaks at Mgi ~ +4.5 and then steadily 
declines toward fainter magnitudes. Our data provide evidence that differential 
reddening may not be ignored in NGC 6791. 

Subject headings: astrometry — open clusters and associations: general — open 
clusters and associations: individual (NGC 6791) — proper motions 



Introduction 



NGC 6791 is an extreme Galactic sta r cluster with an old ag e of ~8 Gyr ( IGrundahl et al. 



20081 ) ■ a high metalhcity [Fe/H]=+0.30 (iBoesgaard et al.l 120091) . and an u n usual orb i t that 



periodically brings it close to the bulge of the Milky Way (IBedin et al.ll2006l ). iKinmanI (119651 ) 
was the first to prov ide an estimate for t he to tal mass of NGC 6791: ~37OOM0 down to 
V =20, confirmed bv iKaluznv fc Udalskil Jl992h . It is widely acknowledged that NGC 6791 
is one of the most massive old open clusters in our Galaxy. 



As indicated by the discovery of several extremely blue subdwarfs (iKaluzny fc Udalski 



19921 ) and the presence of a prominent red clump, morphology of the color-magnitude di- 



agram (CMD) for NGC 6791 is complex. The proposed enhanced mass loss along the red 



Kalirai et al. 


2007 


). Recently, 


Twaroe; et al. 


(2011) 



reported that the CMD in the inner part of NGC 6791 {R < 2') is somewhat different from 
that in its outer part {2' < R < 5'). Specifically, there is a dichotomy near the main-sequence 
turnoff, which the authors interpret as a result of protracted star formation. 

This cluster is critically important to understanding a number of issues, such as stellar 
evolution and population synthesis at high metallicity in the Milky Way and other galax- 
ies, specifically, in ellipticals and in spirals with bulges. Due to its relative proximity and 
its abundance of members, NGC 6791 provides perhaps the best opportunity for detailed 
studies of old and metal-rich populations. Despite a rich array of extant photometric and 
spectroscopic studies, however, the basic knowledge of cluster membership by kinematic 
means has been limited until now t o radial velocities of a few dozen bright stars (e.g.. 



Scott et al. 



1995 



Carraro et al.l 12006) and preliminary prope r motions for a few thousand 



stars dCudworth fc Anthonv-Twaroglll993l : ISedin et al.H2006h . NGC 6791 is located in the 
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field of view (FOV) of NASA's i^Te p/er mission (jBorucki et al.ll2010f) and is one of the targets 
of the Kepler Open Cluster Study flMeibom fc Kepler Teanj|2010l ). It is also one of the clus- 
ters in the WIYN Open Cluster StudjU (jMathieull2000l ) - again, largely owing to its extreme 
properties. 

This Letter presents the highlights of a new and comprehensive proper-motion study of 
NGC 6791, based upon the selection of exquisite deep photographic plates in combination 
with a decade of CCD imaging. The resulting proper motions provide nearly definite cluster 
membership for stars brighter than g' ~22 and, thus, unveil most of the cluster's population 
with unprecedented clarity. 



2. Astrometry and Photometry 



NGC 6791 is a rich but faint star cluster with a main-sequence turnoff at V ~17.5 
(iStetson et al.ll2003l hereafter SBG), thus requiring relatively long exposures even with 3- 
4 m class telescopes. Fortunately, there are deep photographic plates dating back to 1961. 
We selected the best 33 plates (assessed by the depth and sharpness of images) taken with the 
Lick 3 m and the Kitt Peak National Observatory (KPNO) 4 m telescopes. These plates were 
digitised using the Space Telescope Science Institute's GAMMA II multi-channel scanning 
microdensitometer. The second epoch includes 66 frames obtained with the NOAO CCD 
Mosaic Imager at the KPNO 4 m telescope in 1999-2007. The finest 24 CCD mosaic frames 
in g'r' filters (1, 30, and 300 s exposures) were obtained in 2009 November with MegaCam 
at the 3.6 m Canada-France-Hawaii Telescope. These frames served for construction of a 
master catalog consisting of 128,771 objects down to g' ~25 and over ari FOV of 1 deg^. 
We put positions on the system of the UCAC3 catalog (IZacharias et al.l 120101 ) and calcu- 
lated proper motions on this system using the iterative central-plate-overlap algorithm (e.g., 
Jones fc Walkerlll988[ ). The final catalog contains a total of 58,901 objects down to g' ~23.8 
over a smaller area of ~0.8 deg^, defined by a circular FOV (with radius 30') of the KPNO 
photographic plates. The highest precision of our proper motions is 0.08 mas yr^^. Cluster 
membership probabilities are calculated using the frequency functions for cluster and field 
stars, obtained in selected magnitude and spatial ranges. Discrimination between cluster 
members and field stars is excellent down to g' ~22. Details of the astrometric reductions 
and the description of the catalog itself will appear elsewhere (Platais et al., in preparation). 

For photometry, we used only the MegaCam images in g'r' filters. The central part 
of NGC 6791 is fairly crowded, often resulting in blended images. We note that all Mega- 



^This is woes paper 46 of the WIYN Open Cluster Study. 
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Cam images are preprocessed with Elixir which, besides the st andard steps, also norraahzes 
the instruraental photometric zeropoint between the chips (see lMagnier fc Cuillandrell2004j : 



Clem et al.l l2008l ). The short, 1 s exposures obtained in subarcsecond seeing were instru- 



mental in constructing an initial catalog of positions and aperture photometry (MAG_AUTO 
magnitudes from SExtractor). Magnitudes from longer exposures were added only if they 
matched specified tolerances on position and magnitude. In this way it was possible to min- 
imize the deleterious effects of image blending. Nevertheless, image blending and CCD chip 
edge effects may occasionally yield stellar colors that are significantly off. Our instrumental 
photometry was put into the s tandard qV^ syste m by employing the extant secondary photo- 
metric sequence in NGC 6791 (jClem et al.ll2008[ ). The rms error of the calibrating fits is 0.02 
mag. On average, a calibrated g' magnitude near the main-sequence turnoff of NGC 6791 is 
~0.41 mag fainter than that in V . 



New Census of NGC 6791 



One of the most comprehensive photometric surveys of NGC 6791 is that by SBC, who 
also provide extensive lists of bright stars, red giants, and faint blue stars in this cluster — all 
possible or likely members. Now, we are in a position to provide nearly definite star-by-star 
kinematic membership, with a caveat that stars with blended images may not appear in 
our proper motion catalog or, if they do, may be b iased. An exaraple o f such a case is 
the eclipsing binary V20 - a definite cluster member (IGrundahl et al.ll2008l ) but incorrectly 
classed as a field star in our catalog. Its measured proper motion is apparently corrupted 
by the adjacent star (Apos=2") in the amount of ~6cr^. The CMD in Figure 1 shows the 
cleanest and the most complete census of NGC 6791 to date. It contains all stars with 
cluster membership probability >19%, as well as stars with 19% > P^j > 1% but located 
on the main seq uence, all with i n a r adius R < 15' from the cluster center - very close to 
that adopted by iTwarog et al.l ( 120111 ). The cluster membership probabilities for peripheral 
members of NGC 6791 {R > 15') are less reliable because of the overwhelming number of 
field stars and numerical instabilities in finding the cluster population. Hence, in this region 
the selection of cluster members is limited to stars that are on the main sequence and have 
>2%. From the combined CMD we deleted all stars below the main sequence and with 
g' > 19, except for a few very blue stars with P^ > 50%. Likewise, all stars with g' > 16 and 
located well above the main sequence and to the right of the RGB were eliminated. This 
leaves us with a total of 5699 probable cluster members. 
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3.1. Horizontal branch 



The most prominent part of the horizontal branch (HB) — the red clump - — contains 26 



now k inematically-confirmed stars, consistent with the photometric selection by lKaluzny fc Udalski 
(Il992l ). More excit ing is the presence o f extre me horizontal branch (EHB) stars in NGC 6791, 
first recognized by lKaluzny fc Udalski (Il992[). The following EHB stars are kinematic cluste r 
members: B1,...,B7, B9 (names from iKaluzny fc Udalskilll992l : iKaluzny fc Rucinskilll995l ). 



A much fainter star, B8, has a marginal P^=10%; however, its astrometry has large un- 
certainties and therefore we consider it to be a possible cluster member. In addition, star 
B16 (P^=59%) appears to be a cluster member and the following SBG stars are likely clus- 
ter members as well: S 3472 (P^=98%), S13881 (P^=60%). The brightest UV-source, BIO, 
([Landsman et al.lll998[ ) is a definite field star. In Table 1 we provide two new EHB stars in 
NGC 6791, located at P ~ 10'. 

Besides new EHB candidates. Table 1 lists a number of probable cluster members 
that we argue are class ical HB stars. The bluest of these is the well-studied star 2-17 
( jPeterson fc Green! 1 19981 ). which is frequently assigned to the EHB. An intriguing feature of 
the HB here is an upturn at its red end. This could be the effect of duplicity among the red 
clump stars, although that may not be the only interpretation. We stress that proving the 
reality of the HB in NGC 6791 depends critically upon the ability to distinguish HB stars 
from blue stragglers. 

There is one very bright star, S12158 {V =11.96 B — V =1.24), which is formally a prob- 
able cluster member (P^ =91%) but defies classification. Based on its possibly abnormally 
high luminosity, we consider it provisionally to be a field star. 



3.2. Sub-subgiants 



M67 is the first open cluster where the so-called sub-subgiants were discovered ( iMathieu et al. 



20031 ). These apparent binary stars occupy an area of the CMD below the subgiant branch, 
which is not easy to populate with any combination of two normal cluster stars. Owing to 
the low frequency of field stars in this part of the CMD, we are confident that at least five 
likely sub-subgiants are present in NGC 6791. In the notation of SBG, these are: S83, S746, 
S3626, S13753, and S15561. Three of these are also the known variable stars V17, V59, and 



01431_10 listed by lDe Marchi et al.l (120071 ). O nly star yi7 is i n common with the alternative 
list of sub-subgiants (red stragglers) given by iKaluznyl (120031 ). 
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4. Luminosity function and density profile 

The depth and spatial extent of our survey allow the derivation of a reliable luminosity 
function (LF) for the entire cluster down to g' ~ 24. The only omissions are secondaries in 
binary stars and severe blends. We have counted all stars from the cropped CMD (Figure 1) 
in 0.5 mag bins, in two ways: star-by-star and by weighting each star with its membership 
probability (Figure 2). To translate g' into Mg/, an apparent distance modulus of 13.56 
was applied. The peak at Mgi ~ +4.5 and the smooth downward trend to ward fainter 



magni tudes is unmistakable. It is consistent with the appearance of the LF in iKing et al. 



(120051 ). although we cannot replicate t he flattening of their LF at g' > 22 — even less so, the 



flat lu minosity functions provided in iKaluzny fc Udalskil (119921 ) and iKaluzny fc Rucinski 



(Il995l ). A total of ~1000 cluster members per 1 mag bin around the main-sequence turnoff 
{Mg' ~ +4.5) in our present sample is a record high among the known old Galactic open 
clusters. 

We constructed a (number) surface-density proflle for NGC 6791 using the 4830 stars 
in Figure 1 brighter than g' = 22 (where the proper-motion discrimination between cluster 
members and fleld stars is best). This is shown in Figure 3. The flUed circles in this plot show 
the number of stars per unit function of projected cluster-centric radius, deflned by 

adding up stars weighted by their membership probabilities in a series of concentric circular 
annuli. The open circles result from adding up the total numbers of the same stars directly, 
without weighting by P^. 



We fitted a single-mass, isotropic iKing (119661 ) model to the P^-weighted density profile 



in Figure 3, using only the data inside R < 15' to constrain the fit. This implies a rather 
low central concentration for the cluster: c = log(r(/ro) = 0.74 ± 0.05, corresponding to a 
dimensionless central King potential, Wq = 3.4 ± 0.3. The best-fitting King scale radius is 
tq = 4'22 ± 0'30, which, given the low fitted concentration, implies a projected core (half- 
power) radius of Rc = 3'.28 ± 0'12. The associated tidal radius is = 23(1 ± 1(0, and the 
effective (projected "half-mass") radius is Rh = 4(42 ± 0(02 (it is normally the case that 
estimates of Rh are m ore robust than those of other characteristic radii). For a distance of 



4.0 kpc to NGC 6791 (iGrundahl et al.ll2008h . 1' = 1.16 pc; thus, we have that the projected 
core radius of the cluster is Rc — 3.8 pc; the effective radius is P^ — 5.1 pc; and the fitted 
tidal radius is ~ 27 pc. 

Fitting to various other number density profiles defined by both brighter and fainter 
magnitude cuts, and both weighted and not weighted by P^, yields results for Pc, Rh, and rt 
that are generally compatible with the range of values suggested by the numbers just given 
for the g' < 22 case. The same is true for fits to luminosity surface-density profiles defined by 
adding up the luminosities of all the stars in our series of annuli. In exploring these various 
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definitions of the density profile, we saw a clear tendency for the inferred tidal radius (and 
thus the concentration, log(rf/ro)) to increase systematically — along with R^, to a lesser 
extent — as stars with fainter and fainter magnitudes were included to define the observed 
density profile. This is primarily because stars with g' > 19-20 in our sample are somewhat 
depleted, relative to the brighter stars, in the innermost ~2'-3' of NGC 6791. However, it 
is not clear whether this is a physical effect (e.g., due to mass segregation) or an artifact of 
crowding and image blending in the core of the cluster. 

The total luminosity of all stars with g' < 22 in our catalogue is Ltot ~ 5500-6500 i^0,g' 
(assuming D = 4.0 kpc, corrected for Ag/ = 0.55 mag of extinction, and depending on 
whether or not the luminosities are weighted by P/,). Adopting a th eoretical g' stellar mass- 
luminosity relation from the Padova isochrones ( iMarigo et al.l 120081 ) implies a total mass of 
Mtot ~ 5000 Mq for the observed stars. This is a lower limit to the cluster mass because 
w e hav e not a ttempted corrections for stellar binarity or incompleteness. We note that, if 



KingI (119661 ) model with c = 0.74 and Rh = 5 pc accurately describes the internal mass 



profile of NGC 6791, then the projected, one- dimensional velocity dispersion averaged within 
the effective radius should be {a)h — 0.36 {GMiot/ Rh)^^"^ ^ 0.75 km s~^. This corresponds 
to an intrinsic proper-motion dispersion of cr^ ~ 0.04 (D/Akpc)^^ mas yr~^, which is much 
smaller than the formal uncertainties in our proper motions. It is also significantly lower 
than the radial-velo city dispersion of ~ 2 km s ~^ reported for about a dozen red giants by 
Carraro et al.l (120061), suggesting either that the ICarraro et al.l dispersion may be spuriously 
high, or that NGC 6791 may not be in virial equilibrium. Perhaps related to the second 
option, we note that a 5000 Mp^ cluster in a Galactic orbit with a pericenter of ~3 kpc and an 
apocenter of ~ 10 kpc ( jBedin et al.ll2006l ) is expected to have a tidal radius of rt ~ 13 pc ~ 11' 
at pericenter, and rt ~ 28 pc ~ 24' at apocenter — as against a value of rt > 23' suggested by 
our King-model fitting of NGC 6791. Deciding the true dynamical state of this cluster will 
require much more comprehensive and higher-precision proper-motion and radial-velocity 
surveys, to obtain the tightest possible constraints on its orbit and to delineate accurately 
its internal kinematics. 



5. Morphology of the red giant branch 

Even a casual inspection of the cluster's CMD reveals an unusually broad RGB - on the 
order of A{g' — r') ~0.1 mag, while the formal uncertainties in the colors of these stars do not 
exceed ~0.02 mag. There is a variety of potential sources of increased scatter in the CMD, 
such as instrumental effects including the limitations of aperture photometry in crowded 
fields, presence of field stars, differential reddening, binarity of stars, metallicity variations. 
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and possible effects on color /luminosity of RGB stars by variable mass loss suggested by 
Kalirai et all fcoOTh . 



Our proper-motion-vetted CMD offers a better look at the morphology of the RGB, 
which is more sensitive to differential reddening than the main sequence (apart from its 
turnoff region). In addition, our sample covers the largest area ever considered, thus enhanc- 
ing the chances of finding the signs of differential reddening. 

Unfortunately, we cannot take advantage of significantly reduced reddening effects in the 
near infrared. Th ere is a poorly under stood scatter in the J — K color of our RGB stars from 
either 2MASS or ICarney et al.l (120051 ) . which prevents us from using this color as a nearly 
reddening-free reference. Instead, we examined the morphology of the RGB branch in g'r' 
CMD using approximately equal numbers of inner and outer probable cluster members by 
dividing the cluster at R =5'. Then, we draw a ridge line for the inner sample, encompassing 
subgiants and RGB stars; see Figure 4. The majority of RGB stars at i? < 5' are bluer than 
this ridge line but, surprisingly, nearly all RGB stars at i? > 5' are redder than the ridge 
line of the inner sample. The shift in RGB color between the inner and outer sam ples is 
~ 0.05 mag, whereas the age difference of 1 Gyr suggested bv iTwarog et al.l (120111) woul d 
produce only a ~0.01 mag color shift in g' — r' across the RGB region (iGirardi et al.l 12004 ). 
We conjecture that this color shift, at least in part, is due to differential reddening. Among 
the effects listed above, it is the only physical effect known to have an undeniable spatial 
dependency. 

Our data provide evidence that differential reddening may not be ignored across the 
face of NGC 6791. Further advancements in our understanding of such a small effect require 
accurate near-infrared photometry and reliable identification of binaries in the upper part of 
the CMD. 



6. Conclusions 

This study provides a comprehensive and clean list of cluster members in one of the as- 
trophysically most important Galactic star clusters. We hope that our findings in NGC 6791 
and excellent cluster membership will stimulate in-depth future studies of individual cluster 
members. 
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Table 1. New EHB candidates and our selection of HB stars 



ID 


RA (2000) 


Dec (2000) 


9' 


a'- 


-r' 


{%) 


Commenf 


41698 


290? 20761 11 


37? 6103973 


20.291 


-0, 


,140 


46 


EHB candidate 


58783 


290.4511108 


37.7326775 


17.667 


0, 


.032 


56 


EHB candidate 


45106 


290.1308289 


37.6360130 


14.999 


0, 


.508 


27 


S2101 


52460 


290.0726929 


37.6899757 


15.116 


0, 


.845 


42 


8440 


58857 


290.2606201 


37.7331924 


15.186 


0, 


.530 


99 


3-22; RV=-2 var 


60456 


290.2799683 


37.7429771 


15.186 


1, 


.017 


99 


3-27 


64067 


290.1457214 


37.7645378 


15.253 


0, 


.387 


99 


82746 


64589 


290.2497864 


37.7675743 


15.155 


0, 


.278 


99 


2-17; RV=-48 


65895 


290.2140503 


37.7751503 


14.785 


1, 


.075 


99 


NW-20; RV=-54 


69734 


290.2836304 


37.7970619 


14.568 


1, 


.096 


99 


3-33; RV=-20 


69976 


290.1916809 


37.7985649 


14.984 


1, 


.045 


95 


2-45; RV=-53 


82982 


290.3112488 


37.8854256 


15.276 


0, 


.901 


53 


814140 


89107 


290.1376343 


37.9322052 


14.215 


1, 


.101 


41 


82382 


93625 


290.4263306 


37.9681969 


15.165 


0, 


.590 


58 





'^The SBG n umbers are preceded by a letter 8; remaining cross-identifications are from 
iKinmaiil l|l965h : radial velocity (RV) in km is from the literature. 
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Fig. 1. — Color-magnitude and vector-point diagrams of NGC 6791. Upper panel: CMD of 
all probable cluster members, with the exception of some cropped areas. The large dots indi- 
cate our selection of HB stars; the rombs show the proposed sub-subgiants. The arrow shows 



the direc tion and the adopted t otal amount of reddening [E{g' 



0.145, A„> = 0.55, 



following iGrundahl et al.l ( 120081 )]. Lower panel: relative proper motions of probable cluster 
members only. The larger spread of the outer points reflects the relatively low precision of 
proper motions near the limit of the catalog at g' = 23.5. 



- 13 - 




M„ 



Fig. 2. — Luminosity function of NGC 6791. Each bin is 0.5 mag wide. The hatched his- 
togram is obtained by weighting each star, shown in Figure 1, with its membership probabil- 
ity (P^), while the open histogram shows the the actual numbers of stars without weighting 
by Pfj,. The last two bins at absolute magnitude Mgi ~ -|-10 are incomplete. 
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Fig. 3. — Number-density profile of stars in NGC 6791, down to apparent g' = 22. Filled 
circles are the densities obtained by weighting each star with its in concentric circular 
annuli; open circles are the densities obtained by direct add ing up the number of stars. The 
curve is the best fit of a single-mass, isotropic iKingj (119661 ) model to the P^-weighted data 
inside R < 15'. For a distance of 4.0 kpc to the cluster, 1' = 1.16 pc. 
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Fig. 4. — Color-magnitude diagram of NGC 6791 with emphasis on RGB. The left panel 
shows a CMD for the inner cluster sample {R < 5'); the middle panel, for the outer sample 
{R > 5'). For the arrow, see Figure 1. In both panels, the curve indicates the location of inner 
sample's ridge line for subgiants and RGB stars. Note the location of RGB stars relative to 
this ridge line. The right panel elucidates the effect of a reddening E{g' — r') = 0.05 applied 
to the cluster's inner-sample fiducial. This loosely characterizes the degree of differential 
reddening in NGC 6791. 



